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Abstract - This work aims to evaluate the TerraClass 
mapping for the year 2014, in the municipality of 
Paragominas, State of Para. The validation was made hy 
comparing the mapping with the observations found in 
the field. Images of the Satetile Landsat-8, OLI sensor of 
the year 2014, path/row 222/062, 222/063, 223/062 and 
223/063 were used to aid in the field. Using this data it 
was possible to analyze the main representative classes in 
the area, including agriculture, urban area, forest, clean 
pasture, dirty pasture, reforestation, regeneration with 
pasture and secondary vegetation. The secondary 
vegetation presented 2,198.16 km 2 , clean pasture with 
3,332.29 km 2 , agriculture with 896.75 km 2 and the forest 
occupying 54.21% of the total area of Paragominas. The 
overall concordance index was 86%, corroborating the 
reliability of the mapping performed. The average error 
was 6% and the total value of discordance was of 14%. 
Concerning the secondary vegetation, pasture, 
agriculture, urban area and forest classes, they presented 
concordance higher to 50%, while regeneration with 
pasture and reforestation presented greater intensity of 
omission with 40,57% and 76,31% respectively. Inclusion 
errors were less than 40% for the secondary vegetation, 
pasture regeneration, clean pasture and dirty pasture 
classes. The field work was essential to validate and 
analyze the accuracy of the 2014 TerraClass Project for 
the studied region, which becomes important for the 
understanding of the dynamics of land use. 

Keywords - image processing, land use, remote sensing, 
validation 


I. INTRODUCTION 

The spatial occupation and the consequent 
landscape modification have been occurring for decades, 
together with conflicts of interest in various scales. The 
economic development and the expansion of land 
activities have been pointed as one of the main actors in 
the amazon’s deforestation (Feamside, 2006). However, it 
is important to point out that the deforestation does not 
occur as a specific event, but more likely as a process. It 
means that it is not homogeneous in space and time once 
it involves different actors and many productive activities 
as livestock, logging, mining and small and large 
agriculture. Of course, the deforestation is mostly 
detected in the states of Brazil found in the “deforestation 
arch” an area situated in the east and south of the 
Brazilian Legal Amazon (BLA) (Becker, 2005). 

The remote sensing is an indispensable tool in 
landscape or territory scale analysis, once it allows the 
understanding of the degradation processes of a 
determined region (Ferreira et al., 2012). It is used for 
monitoring land use and cover, therefore, it plays a crucial 
role in extracting information from the least accessible 
areas and supporting territory management, due to the 
subsidies that this type of mapping can offer, especially in 
what concern to actions that seek the establishment of 
sustainable practices and the implementation of 
biodiversity conservation policies (Abreu & Coutinho, 
2014). 

The advancement of geotechnologies and remote 
sensing techniques have strongly contributed to society 
with the identification and monitoring of the deforestation 
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(Macedo et al., 2013). State that this set of computational 
techniques has been strongly cooperating for a 
construction of knowledge about spatial and 
socioeconomic patterns. In addition, there are benefits 
from remote sensing application in studies based on the 
integration of field data with land use dynamics and 
vegetation cover (Almeida & Vieira, 2008). 

Studies such as the one realized by Almeida et 
al. (2016) show that with the use of remote sensing, it is 
possible to understand the spatial patterns of land use and 
cover and it contributes to the studies of biodiversity, 
environmental modeling and climate change, essential for 
the creation and monitoring of land use policies. Forest 
formations are monitored by remote sensors in a time- 
scale that ranges from almost daily to annual (Diniz et al., 
2015). 

The conversion of natural landscapes to 
anthropic use or modification of management practices is 
known as land-use change (Foley et al., 2005). Which 
includes selective cutting, commercial planting of trees 
such as eucalyptus and parica, conversion of forest to 
pastures, agricultural production areas, cut-crop 
agriculture and urbanization (Gardner et al., 2009). Which 
has expanded in the last decades, resulting in a 
fragmented landscape with different types of forest cover 
(Laurance et al., 2014). 

The TerraClass Project, executed by the National 
Institute of Space Research (INPE), in a partnership 
between the Amazon Regional Center, Embrapa Oriental 
Amazon and Embrapa livestock and Agriculture 
Informatics has as main objective qualifying the 
deforestation of the Legal Amazon, based on the 
deforestation areas mapped and published by PRODES 
Project and satellite images (Inpe, 2016), playing an 
important role in determining which land uses are 
replacing the Amazon forest (Almeida et al., 2016). In 
addition, this project is motivated by the concern of the 
scientific world and society in general with the threats to 
the greater biodiversity area of the planet, considering the 
role of Amazon forest in the context of global climate 
change (Inpe, 2016). 

The interest in monitoring forest resources has 
been increasing by the past decades and so the demand 
for mapping plant cover at regional and global scales 
(Shimabukuro, 2000). Thus, validating orbital remote 
sensing data through the verification of truth land points 
(Congalton e Green, 1999), has become difficult and 
sometimes impractical in some parts of the (Hess et al., 
2002 ). 

The Para State has the largest deforested area of 
the BLA, with 262.088 km 2 , equivalent to 33.72% of the 
Amazon’s total deforested area until 2016 (Inpe, 2017). 

This scenario is due to the installation of agromineral 
project in its territory, as the “Grande Carajas” Project for 
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iron ore exploitation, the implementation of transport and 
energy infrastructures and the advancement of the 
agricultural frontier (Kohlhepp, 2002). 

The accuracy assessment of the land use and cover 
maps has been well studied by researchers in remote 
sensing, however, the field measurements of large 
territorial mappings present obstacles because of the 
difficulty of access, operational and logistic costs, as well 
as the high planning and execution demanding time 
(Adami et al., 2012). 

According to Almeida et al. (2016), the first 
TerraClass mapping had a global accuracy of 76%, when 
tested to the states of Para and Mato Grosso, attending the 
accuracy’s need of the mapping and allowing statistic 
inferences (McRoberts, 2011). 

In this sense, this papes has the objective of assessing 
the TerraClass mapping for the year of 2014, in a local 
scale, which will be applied to the municipality of 
Paragominas, in the state of Para. 

II. WORK METHODOLOGY 

The study was realized in the municipality of 
Paragominas, situated at the margins of the Belem- 
Brasilia Highway (BR-010), with connections with the 
PA-125 e PA-256 state highways. The population is 
estimated in 97.819 inhabitants, distributed in an area of 
19.342.254 km? (IBGE, 2010). The municipality is part of 
the southeast region of Para State and the Paragominas 
micro region. It presents limits with the municipalities of 
Ipixuna do Para e Ulianopolis and its municipal seat is 
limited by the following geographic coordinates: 
03°00’00” S e 47°21’30” W (Fig.l). 



The utilized images are from the Lands at-8 OLI 
(Operational Land Imager) satellite, which are available 
free of charge from the USGS (United States Geological 
Survey), with spatial resolution of 30 m. The images were 
from 2014, of the points 223 and 222 and orbits 62 and 
63, dated from June, October and November. 

The images were selected based on the lower 
presence of clouds and then submitted to geo referencing 
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and contrast enhancement treatment. At last, it was 
realized the procedures of equalization, to obtain better 
identification, and grouping of the land use and cover 
features (Fig. 2). 



Fig.2: Path/row of the 2014’s Landsat-8 images 


The municipality of Paragominas was chosen 
due to the fieldwork that was executed in 2014, through 
Sustainable Landscapes Project, which involved a land 
route and an overflight to register images of 
georeferenced points, generating pictures that based the 
mapping validation. 

The validation of the land use and coverage map 
for the year 2014 was based on fied data, collected with 
GPS navigation device for land and air. In total, 321 
points were collected in the field, which were used to 
calculate global accuracy (Fig. 3). 



Fig.3: Fieldwork visited points in the municipality of 
Paragominas, PA 

The accuracy assessment is an important process 
in the remote sensing data analysis because it results in 
statistic inferences to the reliability of the generated data 
(McRoberts, 2011). 

The overall accuracy is one of the simplest 
measures, calculated by the total sum of correctly 
classified pixels divided by the total number of pixels in 
the confusion matrix. In addition, individual category 
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accuracy is calculated through the accuracy of the 
producer, indicating the likelihood of a reference pixel be 
correctly classified, and the user's accuracy, which 
computes the probability of a mapped pixel be compatible 
with the field class (Congalton, 1991). 

To calculate the Global Accuracy or Global 
Precision (GA), User Accuracy (UA) and Producer 
Accuracy (PA), the georeferenced pictures taken in the 
fieldwork were analyzed together with the TerraClass 
2014’s mapping for the study area. For the calculation of 
accuracy, a contingency table was built (Pontius & 
SantaCruz, 2014) (Table 1). 


Table.1: Example of a contingency table containing the 
functions for calculating the GA, UA e PA. 




Reference Classes 

Total 

Producer’s 



Classi 

Class... 

Classc 

Accuracy 

Mapping 

Classes 

Class i 

Pn 

Pi... 

Pic 

Pi+ 

PAi=Pii/Pi + 

Class... 

P...i 

P 

P...C 

P...+ 

PA = P IP + 

Classc 

Pci 

Pc... 

Pec 

Pc+ 

PA c =Pcc /Pc+ 

Total 

P + 1 

p + 

Pc + 

n 

User’s Accuracy 

UAi=Pn/P +1 

PA... = 

P /P + 

PAc=Pcc /Pc+ 

“-I? 


The values allocated on the main diagonal of Table 1 
represent the correctly classified elements, whereas the 
elements outside the main diagonal represent errors of 
omission and inclusion (Stehman and Foody 2009). 

III. RESULTS AND DISCUSSIONS 

The technique of visual interpretation of the spectral 
targets of a given satellite image is part of data generated 
through Remote Sensing. In this sense, the acquisition of 
points of terrestrial truth has become an important tool for 
validating the data generated by interpreters through 
orbital images (Espifito-Santo, 2005). 

The TerraClass Project maps the following 
thematic classes: Annual agriculture, Mosaic occupations, 
Clean pasture, Dirty grass, Regeneration with pasture, 
Pasture with exposed soil, Secondary vegetation, 
Reforestation, Urban area, Others, Mining and Area not 
observed (Adami et al., 2015). With the exception of the 
annual Agriculture class, all other classes are mapped 
through photointerpretation. The annual Agriculture class 
is mapped from an automatic method based on the 
spectral-temporal behavior of the Normalized Difference 
Vegetation Index (NDVI), obtained by the MODIS sensor 
(Roerink et al., 2000). 

The validation for the municipality of 
Paragominas was carried out with the following thematic 
classes of mapping of the TerraClass Project, year of 
2014: agriculture, urban area, forest, clean pasture, dirty 
pasture,reforestation, pasture regeneration and secondary 
vegetation, representative in the municipality. 

Based on the TerraClass 2014 mapping, the 
forest occupies 54.21% of the municipality, with 
10,483.37 km 2 . The secondary vegetation presented 
2,198.16 km 2 , the clean pasture with 3,332.39 km 2 , the 
annual agriculture with 896.75 km? (Fig. 4). 
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Fig. 4: Use and land cover map year 2014 of the 
municipality ofParagominas , PA. 


In all, these three classes correspond to 
approximately 85% of the area of the municipality. 
Another interesting factor to be observed is that the 
predominance of forest is inversely associated with the 
distance of the roads that cut the municipality. 

From the satellite image and photographs, it was 
possible to represent some validation points that can also 
be used for training in photointerpretation. 

In point 1 it was possible to observe and validate 
in the field the class of urban area, with characteristics of 
streets, squares, houses positioned in a very closely and 
with regular spatial distribution (Fig. 5). 


represented with soybean planting, in the magenta color, 
smooth texture and regular polygons characterize the 
feature (Fig. 6). 



Fig. 6: Class of Annual agriculture. 


It was observed in section 3 the class clean 
pasture validated in the field according to the 
representation of well managed pasture area, with low 
infestation of herbaceous and shrubby weeds (Fig. 7). 



Fig. 5: Class of Urban area. 



In section 2 the validation was verified by the 
observation in the field of the annual agriculture class 
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For point 4, the regenerativo pasture class was 
confirmed with the field truth by the presence of invasive 
species at various levels of development (Fig. 8). 



Fig. 8: Regeneration pasture class. 


For point 5 the validation of the Reforestation 
class, represented by areas with regular polygons, division 
of plots and homogeneous texture (Fig. 9). 



Fig. 9: Reforestation class. 


At point 6 the field truth was verified by the 
Forest class represented by the forest management reserve 


area, in the satellite image with rough texture and dark 
green color characteristics (Fig. 10). 



Fig. 10: Forest class. 


Point 7 represents the dirty pasture class with 
presence of many herbaceous and some shrub species 
(Fig. 11). 



Fig.11: Shrubby pasture class. 


For the representation of point 8 the secondary 
vegetation class was validated proving the aspect of 
natural regeneration of native shrub and tree vegetation 
(Fig. 12). 
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Fig. 12: Secondary vegetation class. 


The classification accuracy is expressed using 
the error matrix also called the confusion or contingency 
matrix, derived from the classification and reference data 
set as the starting point of the precision measurements 
(Story and Congalton, 1986), which shows the cross 
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tabulation of the mapped land use and cover versus the 
actual one assessed in the field (Foody, 2002). The 
contingency matrix lists land use and cover reference 
values in the columns and the classified data in the rows, 
being the main diagonal of the matrix the correctly 
classified data (Banko, 1998). 

Based on the contingency matrix, global 
accuracies, omission and inclusion errors per class and 
indices of disagreement (displacement, quantity and 
change) were calculated (Pontius and Santa Cruz, 2014). 
Studies performed in Cerrado areas show high reliability 
with a general agreement index of 80.2% (Maurano and 
Adami, 2017). 

The validation process is a necessary step in a research 
that aims to analyze in a statistical way the accuracy of 
the classified images. Thus, table 2 shows the results of 
the contingence matrix of the classification for the 
municipality of Paragominas. 


Tahle.2: Comparative Results between mapped data and field data 


TerraClass 

mapping 

class 

Agriculture 

Urban 

area 

Forest 

Clean 

pasture 

Shrubby 

pasture 

Reforestation 

Regeneration 
with pasture 

Secondary 

vegetation 

Total 

Agriculture 

50 



1 


3 



54 

Urban area 


7 







7 

Forest 



110 

3 



1 

3 

117 

Clean 

pasture 

9 



49 

2 

12 

3 


75 

Shrubby 

pasture 




1 

7 

1 

2 

1 

12 

Reforestation 






7 



7 

Regeneration 
with pasture 







7 

6 

13 

Secondary 

vegetation 



4 

3 


4 


25 

36 

Total 

59 

7 

114 

57 

9 

27 

13 

35 

321 


The general agreement index between mapping 
and the validation was 86%, demonstrating high 
reliability of the mapping performed by TerraClass, based 
on the result of the contingency matrix. 

According to Pontius & Santa Cruz (2014), the 
mean values of agreement and disagreement are divided 


into quantity, displacement, and change. The quantity 
errors refer to the absolute difference between omission 
and inclusion errors, presenting an indicative of the mean 
error of the mapped area. The displacement refers to the 
unpaired allocation difference between the analyzed 
classes in the contingency matrix and, at last, the change 
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refers to the error caused by paired exchanges between 
classes in the contingency matrix (Foody, 2009). 

Thus, the mean error presented for the area was 
6%, displacement of 7% and the change of 1%. When 
these disagreements are summed, a total error of 14% is 
obtained, which is considered low due to the complexity 
of the mapped classes and size of the study area (Fig. 13). 
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” CONCLUSION 

The presentation of the contingency matrix is 
fundamental, since it is possible to visualize the 
confusions that occurred between classes. 

The global accuracy may be sufficient, depending on the 
purpose of the map. However, if there are specific classes 
of interest or more than others, individual classes’ 
accuracy may be relevant. 

The global accuracy of the TerraClass mapping for the 
municipality of Paragominas was 86%. Therefore, future 
research on this project should consider improvements for 
the mapping of regeneration with pasture and 
reforestation classes. 

The validation, as part of the mapping process, estimates 
the data accuracy of the land cover use dynamics, giving 
support to the information consolidation. 
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